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1. Introduction 



A precise calculation for neutral B mixing matrix elements has been an urgent task since the 
recent measurement of Am., [jl|]. Such a calculation is crucial in obtaining stringent constraints on 
the unitarity triangle of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and searching for new 
physics. On the other hand, the D" mixing system is a good channel to search for new physics 

because the Standard Model contribution is strongly suppressed. In the Standard Model, the 
short distance contribution to the mass differences of these heavy neutral meson mixing systems is 
predominantly determined by the matrix elements of a single set of four quark operators: 



'l,aa 



where /i is a heavy quark field (either a ft or a c quark), qa is a light-quark field with flavour a {a is 
not summed over), and a and j3 are colour indices. Models containing flavour-changing currents 
other than the V — A form (arising in supersymmetric extensions of the Standard Model and other 
scenarios) usually result in mass differences that additionally depend on matrix elements of the 
four-quark operators ||3|| 

ff2,aa=h''{l-Y5)q^ahH^-Ys)ql 

^Xaa = /^"(l-75kf /l^l-75k«", (1-2) 

G,^aa=h"{\-rs)qrh^l + r5)q^, 



Generically we can represent these operators as 

Giaa=hriqhr2q, (1.3) 

for the appropriate choice of spin and colour matrices, Fi 2. In lattice calculations, it is convenient 
to perform a Fierz transformation which renders linear combinations of the operators in Eq. ( pr2| ) 
into products of colour-singlet currents. Nevertheless, we choose to work in the basis of Eq. ( |1.2D . 

In this article, we present a study for the light-quark mass and spatial volume dependence in 
matrix elements relevant to the mass differences in the neutral B and D meson mixing systems, 
using (partially quenched) heavy meson chiral perturbation theory ^ So far, only one exploratory 
numerical calculation in quenched lattice QCD has been carried out for the full set of these matrix 
elements [||]. However, with experimental progress in CKM physics, an accurate determination of 
these matrix elements, which involves reliable chiral extrapolations, will become necessary in the 
foreseeable future. In this work, we show that the chiral extrapolation for matrix elements of Gi ^a 
is considerably less complicated than that for matrix elements of the operators in Eq. ( |1.2| ). This 
is due to the fact Gi^aa preserves heavy quark spin symmetry, as explained in Sections |3| and ^. 
Details of the formulae for the chiral extrapolations of these matrix elements have been presented 
in Ref. [§. 

' Some of the matrix elements studied in tliis work are also relevant to the width difference in the B mixing sys- 
tem [§]. 
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2. Heavy meson chiral perturbation in finite volume 

Heavy meson chiral perturbation theory (HM;^PT) was formulated in Refs. [Q, ||, ^, and gen- 
eralised to the quenched and partially-quenched versions in Refs. |1C, 11]. Finite spatial volume 
effects in this effective theory have also been investigated in Ref. [12]. Here we only sketch the 
ingredients necessary to study neutral heavy-light meson mixing systems. Details of this effective 
theory can be found in the above references. 

In HM;^PT, the heavy-light meson field appears in the covariant form 

H^^'='-±^{p*S'f-P^Fh), (2.1) 

where Pa^^ and P*i[^^ annihilate pseudoscalar and vector mesons containing a heavy quark Q and 
a light anti-quark of flavour a. In the heavy particle formalism, such mesons have momentum 
=Mpv^ with l^'^l <^Mp and is the velocity of the particle. Under a heavy quark spin 
SU (2) transformation S and a generic light-flavour transformation U [i.e., U G SU (3) for full QCD 
and U G 5C/(6|3) for partially-quenched QCD (PQQCD)], 

hP^SH^^^uI. (2.2) 

The conjugate field, which creates heavy-light mesons containing a heavy quark Q and a light 
anti-quark of flavour a, is defined as 

HiQ) = fH^Q^y, = (p;(e)t^. +p(e)t^^) [+1^ (2.3) 

which transforms under S and U as 

Hi^^ Ua,Hfs\ (2.4) 

These heavy-light meson fields couple to the Goldstone meson fields 

^ = e'"^/^ = ^/l, (2.5) 

where Z is the ordinary nonlinear Goldstone field. The field transforms as 

^^Ui^^U"^ = U^UI (2.6) 

where t/L(R) element of the left-handed (right-handed) SU{3) and SU {6\3) groups for QCD 
and PQQCD respectively. 

The mesons containing a heavy anti-quark Q and a light quark of flavour a can be included in 
the theory by applying the charge conjugation operation to the above heavy-light meson field //i^^. 
The field that annihilates such mesons is 

//F = (<f7'^-Pi^V5)^, (2.7) 
which transforms under S and U as 

hP ^UabH^^^S''. (2.8) 
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The main effects of heavy quark symmetry breaking in HM;^PT Lagrangian is the splitting 



between vector and pseudoscalar heavy-Ught meson masses 1 13] 



— tro \^Ma O^yMa 



(2.9) 



where A2 is a low energy constant (LEC), and trn means trace in the spinor indices. 



3. Four-fermion operators 



Under a chiral transformation, the four-quark operators in Eq. (1.2) fall into two categories: 



5ll = h Fll qt h Fll qt, 
?LR = h qi h r[^' qR, 



(3.1) 



where 



Ql,r 



1±75 



-q. 



(3.2) 



Operators ^i^aa, ^2,aa and ^3^aa are of the first type and transform in the symmetric (6l,1r) 
representation built from the direct product (3l, 1r) <8> (3l, 1r) = (6l, 1r) © (3l, 1r) under chiral 
rotations while i&4_aa and ^s.aa are of the second type and transform in the (3l,3r) representa- 
tion. Here we refer to the SU(3) flavour transformation properties, leaving the partially quenched 



extension to the following subsection. Note that the colour indices in Eq. ( |1.2| ) are relevant to short- 



distance physics, and hence play no role in the chiral properties of these operators. Treating r^^, 
and r[^^ as spurions transforming as 



r(i) 

^ LR ' 

r(2) 

^ LR 



S Tll U, 



t 



^ LR "^fi ' 



(3.3) 



the operators in Eq. (3.1) remain invariant under heavy-quark spin and chiral rotations. The boson- 



isation of the operators in Eqs. (1.1) and (1.2) can be performed using these spurion transformation 
properties. The procedure is explained in details in Ref. [^]. It leads to the following set of opera- 
tors involving the individual heavy meson fields: 



^HM;t;PT 

I. cm 

5,HM;t:PT 



'2(3].aa - P2{3) [^P 



yHM;i;PT 

4{5),aa 



^4(5) {^P"^' 



(3.4) 



where the pi, j3/,j3,-, and j3/ are LECs. 

It is important to note that in the above equation, the operator ^f^J^^'^ behaves somewhat dif- 
ferently from the other operators as only a single LEC, jSi, occurs. This greatly simplifies any chiral 
extrapolation of corresponding lattice data for neutral heavy-light meson mixing in the Standard 
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Model, as discussed in the next section. We stress that this simplification is not obvious from the 
operator structure and is particular to the V — A form of the Standard Model currents. In general, 
one would expect that operators for pseudoscalar and vector meson mixing processes are accom- 
panied by different LECs. This is the case for all the non-Standard-Model operators, as shown in 
Eq. (p 



To understand the origin of the above simplification in the Standard Model operator i^™^^^, 
we turn to heavy quark effective theory (HQET). In this effective theory, the operators that produce 



the same matrix elements as those in Eq.(|1.3|) are [ ]14| ] 



aHQET 



(3.5) 



where Fi 2 are the appropriate Dirac and colour structures from Eq. (1.2). Here, Q and Q denote 
fields annihilating a heavy quark and heavy anti-quark, respectively (these fields do not create the 
corresponding anti-particles). Additional terms in HQET which create two heavy quarks or anni- 
hilate two heavy anti-quarks will not contribute to neutral heavy-meson mixing and are ignored. 
The standard model operator in HQET, i^f2f^' satisfies the relation 



{4><r}i^)={4<r}i^) = o, 



(3.6) 



where \P) is pseudoscalar heavy-hght meson state, and \P) is the state of its anti-particle. The 
operator 

Sl = e'i\Q^OijQ-QaijQ\ (3.7) 



is the heavy quark spin operator [ ]15P that changes the spin of the heavy-light meson state by 
one. Therefore, Eq. ( |3^ ) implies that the mixing matrix elements for vector and pseudoscalar 
heavy-light mesons are equal and opposite in the heavy-quark limit. This symmetry is reflected in 
HM;t:PT, leading to the result for in Eq. (gj)- 

For the non-Standard-Model operators, it is straightforward to show that 



and 



r,3 ^HQET 



r,3 ^HQET 



(3.8) 



(3.9) 



where / = 2,3,4,5. This means that the pseudoscalar and vector meson mixing processes via these 
operators are not proportional to each other, hence the appearance of the terms accompanied by 
3 4 5 and jS^j in Eq. 

We end this section by noting that equations of motion for the heavy quark [Q] result in ^3,00 = 
— ^i,aa/2 — &2,aa, and Can be used to relate some of the LECs in Eq. (3.4). 



4. Neutral meson mixing matrix elements 

Calculations at NLO in the chiral expansion require the evaluation of the one-loop diagrams 
shown in Figure [l[ We perform these calculations both at infinite volume and in a cubic spatial 
box of dimensions (the time extent is assumed to be infinite). In this section we summarise the 
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/ 



(a) 



(b) 



(c) 



Figure 1: Diagrams contributing to the matrix elements of four-quark operators at NLO in the chiral ex- 
pansion. Solid, double and dashed lines correspond to propagators of pseudoscalar and vector heavy-light 
mesons, and Goldstone mesons, respectively. The crossed circle denotes the four-quark operator and dia- 
gram (a) is the wave-function renormahsation. 

results, relegating details of the calculations to Ref. [^. Although we present results specifically in 
the B-meson systems, note that they are also applicable to Z)-meson systems, under the assumption 
that the charm-quark mass is large enough compared to Aqcd- 

For the standard model operator we find the following matrix elements 



+ ) + analytic terms , 



(4.1) 



'l,vs 



j8i( 1 + 5-/'' + 



^5? + 



+ ^i^' ) + analytic terms . 



The "analytic terms" here include Goldstone meson mass squared terms, a term ~ as{M},)/An 
(arising from mixing) and a term ^ Aqcd/^/> The wave-function contributions, Wm (diagram 

(a) in Figure |l|), and tadpole- and sunset-type operator renormalisations, =3^/'' and ' (diagrams 

(b) and (c) in Figure respectively) are non-analytic functions of the light quark mass and lattice 
volume and are defined in Ref. [Q]. 

For the operators that contribute to the B-meson mixing processes beyond the Standard Model, 
we obtain: 



,0| 



'2{-i),ss 



\B" 



0\ 



^2(3) 
^2(3) 



, «/ ^(2(3)) 



(4.2) 



0\ 



{B \^4{5),dd 
' dO I I dO 



^4(5) + ^4(5) 



1 + =^. 



(4(5)) 



+ ■ 



{B, I ^4(5),i.5 \B^ ) = /34(5) + ^4(5) 



1 + + 



2 



+ 



+ 



^4(5) + ^4(5) 

^4(5) + ^4(5) 

(0 



(4(5)) 



(4(5)) 



+ ... 



+ ... 



where ". . ." represents the "analytic terms" as in Eq. (4J )■ The terms ~ arising from the sunset 
diagrams [Fig. [l](c)] involve the neutral heavy-light vector meson mixing amplitudes. As discussed 
in the preceding section, it is only in the case of ^i^qq that these amplitudes are related to those of 
the pseudoscalar heavy-light mesons. For / = 2,3,4,5 these terms are consequently accompanied 
by different LECs. 
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The above results lead to the fact that the chiral extrapolations for B mixing matrix elements 
in and beyond the Standard Model have different features. Generically, the chiral expansion for 
{BlWi^aa\Bl.) takes the form 



{bI^ I &ua\Bl^) 71 (1 + ^) + analytic terms, (4.3) 



where 71 is the leading-order LEG, L denotes the non-analytic one-loop contributions (chiral log- 
arithms), and the analytic terms are from the next-to-leading-order counter-terms in the chiral ex- 



pansion. However, for the operators in Eq. (|1.2|), the chiral expansion has the generic feature: 

(5° )) 7 (1 +L) + iU + analytic terms, (4.4) 

where / = 2,3,4,5, 7 and are unknown leading-order LECs, and L and L' are different one- 
loop chiral logarithms. The appearance of the second non-analytic term complicates the chiral 
extrapolation in Eq. ( p!4| ) because an additional unknown parameter must be determined. 

We end this article by noting that the effects of scalar resonances in these matrix elements have 
been studied in Ref. [|l6|] and found to be negligible. 
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